A Kohn-Sham system at zero temperature 
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1 Introduction 



Hohenberg and Kohn have shown in [7j that the ground state of an iV-body 
quantum system at zero temperature is completely determined by the particle 
density Nowadays that paper is considered as the starting point of the so- 
called density functional theory. The main advantage of this approach is 
that the description of an iV-body quantum problem can be reduced to an 
effective one-body system. A shortcoming of [7] is that only the existence 
of such a functional, depending on the particle density, was shown, but the 
functional was not given explicitely. In |14| Kohn and Sham have indicated 
such functionals which ^-electrons look like 

-2 N 



n=l J 



; ° k=l 

where M is the number of positive ions, Zf, their atomic number, the 
positions of ions, q is the magnitude of the elementary charge, e is the 
vacuum permittivity and e xc [u] is the so-called exchange correlation energy 
density. The particle density u is given by the expression 

N 

M (r):=2^|^ n (r)| 2 , (1) 

n=l 

where 2 counts for the spin degeneracy of the particles and ip n are eigenfunc- 
tions satisfying the Kohn-Sham equation 

- /- J2 + q* [ - ^1 dh> + V xc [u](r)) i> n = E n ^ n . 

2m 47re ^|r-R fc | J 4vre |r - r'| J 

By 14 c [w] the so-called exchange correlation potential is denoted, which is 
given by 

V«[ti](r) : = ^MMi. 

ou 

The potential 

q 2 M Z 

Va(T) ~- = -^Xv^Y 

fe=l 

which is determined by the positive ions, can be regarded as a given external 
potential. The potential 

/ \ f U ( r ') j3 I 
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is nothing else as the solution of the Poisson equation 



Mr) = 



(2) 



So we end up with a Schrodinger operator of the form 




(3) 



with the effective Kohn-Sham potential 

V := V + V xc [u] - qf, 



(4) 



where if obeys the Poisson equation (j2j). The model is very flexible and 
widely applicable because the exchange correlation term can be well adapted 
to a great variety of problems. This is one of the reasons why the approach 
was very successful in the last forty years, and Kohn was awarded the Nobel 
prize in 1998 for that idea. 

We note that the Kohn-Sham system is quite similar to the so-called 
Hartree-Fock approximation in iV-body quantum systems [TSJIlH]- The main 
difference is that the exchange correlation term for the Hartree-Fock system is 
not local. However, performing the low density limit one obtains the Hartree- 
Fock-Slater approximation [231 HI 12]- in this case the exchange correlation 
potential is of the form V^ c [m](z) = —C a \u(x)\ a , a = 1/3 in three dimensions 
(3D). For a = 2/3 one gets another interesting approximation which is called 
the Thomas-Fermi correction. Usually, models of that type are summarized 
as Schrodinger- Poisson-X a systems, see [TJ, which fit, of course, into the 
class of Kohn-Sham systems. In the following we do not restrict ourself 
to Schrodinger-Poisson-X a systems, but consider a larger class of local and 
non-local exchange correlation terms including the X°-models. Generally, 
we assume that the exchange correlation is a non-linear mapping acting from 
the set of densities into the set of potentials obeying a certain continuity 
condition. 

Note that all these considerations are made at zero temperature. An 
extension of the Hohenberg-Kohn approach to temperatures above zero was 
proposed by Mermin in [IT]. He showed that the expression for the particle 
density then modifies to 




n=l 



such that 
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where (3 := 1/kT and fi is the so called chemical potential. 

There are many papers on the numerics of the Kohn-Sham system, but 
very few on its mathematical analysis. In case of nonzero temperature and 
bounded domains the system was analyzed in [9] and [TO], where existence 
and a priori estimates were shown. In [20] the Schrodinger-Poisson-Slater 
system was investigated for a periodic external potential V . The time dy- 
namics of the Schrodinger-Poisson-Slater system was considered in |22j. 

In the following we are going to investigate the zero and non-zero tem- 
perature Kohn-Sham system with a general exchange correlation potential 
for a planar semiconductor nanostructure. The system reduces essentially to 
an effective one-dimensional system. Since for one-dimensional systems the 
eigenvalues are simple, one avoids in this way the occupation problem for 
the last eigenvalue at zero temperature, if it is degenerated. We show the 
existence of solutions for such systems at non-zero and zero temperature. In 
particular, we prove that the solution is unique, if the exchange correlation 
potential is absent. In the zero temperature case this proof is based on an 
extension of the monotonicity for the negative particle density operator to 
non-smooth distribution functions like Fermi-Dirac distribution function at 
zero temperature, see [Sj. Finally, we prove that the non-zero temperature 
solutions of the Kohn-Sham system converge to those for zero temperature 
as the temperature goes to zero. 

The outline of the paper is as follows: In Section 2 we derive an expres- 
sion of the effective one-dimensional particle density for a planar semicon- 
ductor nanostructure. In Section 3 we introduce the mathematical setup of 
an one-dimensional Schrodinger-Poisson system and make it mathematically 
rigorous. Section 4 is devoted to the existence of solutions. In Section 5 we 
prove the uniqueness of solutions, if the correction term is absent. Finally, 
in Section 6 we show the convergence of non-zero temperature solutions to 
zero-temperature ones as the temperature goes to zero. 

Notation: In this paper the system is considered on the domain Q := 
]0, 1[. For this reason we omit for all functional spaces the explicit indi- 
cation of this interval; e.g. write L l instead of -^ 1 (]0, 1[) and so on. We set 
Wq' 2 : = {/ G W 1 ' 2 : f(0) = f(l) = 0}. The space of antilinear forms on W 1,2 
is denoted by W~ 1,2 . For Banach spaces X and Z, we denote by B(X; Z) the 
space of all linear, continuous operators from X into Z. If X = Z we write 
B(X). Because of the numerous use of X = L 2 we introduce the abbreviation 
II ' II = II ' \\b(l*)- 
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2 Particle density for planar nanostructures 



We consider a planar semiconductor nanostructure; that is, there is a se- 
quence of layers of different materials along the x-direction (i.e. a sequence 
of quantum wells and barriers) embedded between two thick layers of isolator. 
Then the wave functions of a particle (electron or hole) are given by 



e ik x r x 



W r ) = (7) 
and the total energy of the particle is 

E = ¥± + X l , (8) 
2m± 

where r_|_ = (y, z) represents the transversal coordinates, kj_ = (k y , k z ) the 
transversal wave number and m± the effective mass in the transversal di- 
rection. The wave functions along the x-direction and their corresponding 
energies correspond to the one-dimensional Schrodinger operator in the ef- 
fective mass-approximation (Ben-Daniel-Duke form) 

H v :=--—[m- 1 —)+V, (9) 



2 dx \ dx / 

where m — m(x) is the position dependent effective mass and V is an effec- 
tive Kohn-Sham potential to be specified later on. The embedding isolator 
layers impose homogeneous Dirichlet boundary conditions for the wave func- 
tion along the x-direction, providing a discrete spectrum of energies, A;, and 
defining in such a way a closed system in the x-direction. 

Quantum mechanically, the particle density is given by a sum over all 
states of their localization probability multiplied by their occupation proba- 
bility. For fermions (electrons and holes) the occupation probability is given 
by the Fermi-Dirac function. 

At zero temperature, all states up to the Fermi energy E F are equally oc- 
cupied, with probability 1 and above Ep all states are empty, i.e. occupation 
probability 0. Thus, for zero temperature, the particle density is calculated 
quantum mechanically as 



u(r) = 2 / dk ± J2\^4 
J i 



E<E F 



where 2 counts for the spin degeneracy of the particles. Using the expression 
(J7J) of the wave functions for planar structures, one has only an x— dependent 
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particle density 

where the sum runs up to the last occupied level, i.e. Xn f < Ep and the 
integral is taken up to a maximum value of the transversal wave number, 

k± F = \J^jp~ (Ep — Xi) - depending on I. The integral over <ikj_ can be 
performed, and one obtains the particle density at zero temperature for an 
effective one- dimensional system [21] as 



N F 

u(x) = 2^- 2 J2\M^)\ 2 (Ef-\i). (10) 
1=1 

At temperature T different from zero the particle density is given by 
w(r) = 2 J dk x J2\^( r )\ 2 ME,l^) 



2 / dk ± V|vI/ kx>/ (r)| 2 /FD(^/i), 
Jo i=i 



where fpD{E,ij) is the Fermi-Dirac distribution function 

f FD (E,fi) = — , 

1 + e kT 

k is the Boltzmann constant and \x is the chemical potential. Inserting ((7|) 
and using (jHJ) one obtains an x— dependent particle density 

This corresponds to the general form (f!6l) used later on. Also in this case 
one can perform the integral over dh±, obtaining [31 125] 

oo 

u(x) = 2^- 2 kT \Mx)\ 2 In (l + e^) , (11) 
71 i=i 

which provides the expression for fp(x) used in Lemma [6.51 below. Carrying 
out the limit T — >• in (fTT]) one obtains the same expression for the particle 
density as in the zero temperature limit ffTUj) . because 

lim kT In (l + =(p- A,) G(p - A/) 



(G being the Heaviside function) and it is known that the chemical potential 
for zero temperature equals the Fermi energy, lim fi = Ep. 

3 Kohn-Sham systems 

3.1 Kohn-Sham systems in one dimension 

The Kohn-Sham system is a system of equations governing the electrostatic 
potential p> and the particle density u under consideration. Let us consider a 
system of electrons. The electrostatic potential and the particle density have 
to obey Poisson's equation 

-1^ = °-^ (12 > 

in the device domain Q = (0, 1) where q is the magnitude of the elementary 
charge, and e = e(x) denotes the dielectric permittivity. The right-hand side 
of (TT2l) is a charge distribution D of ionized dopants and the particle density u 
which is defined below, see ( Fl6l) . One has to supplement the Poisson equation 
(112j) by boundary conditions. Usually one chooses inhomogeneous Dirichlet 
boundary conditions 

<p(0) = ^ 6l and <p(l) = (pi e M, (13) 

which model Ohmic contacts. A straightforward calculation shows that inho- 
mogeneous boundary conditions can be transformed into homogenous bound- 
ary conditions. Indeed, introducing the function p> : [0, 1] i — > M, 

and setting (f> := <p(x) — (p, x G [0, 1], one gets that (p satisfies the Poisson 
equation 

d d , 

- t x £ tJ = D ~ qu (14) 

and obeys the homogeneous Dirichlet boundary conditions 

0(0) = and 0(1) = 0. (15) 

This gives rise to the Poisson operator V := supplemented by ho- 

mogeneous Dirichlet boundary conditions. 
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The particle density u is computed by the quantum mechanical expression 

oo 

u(V)(x) = 2 £ /(A,(V) - |^(V)(ar)| 2 , x G (0, 1), (16) 

z=i 

where 2 counts for the spin degeneracy, / is a distribution function (to be 
specified later on), A; = \i(V) are the eigenvalues and ipi = tpi(V) are the 
corresponding L 2 -normalized orthogonal eigenfunctions of the Schrodinger 
operator H v , cf. The chemical potential /i/(V) is determined by the 

condition 

oo 

2^2f(MV)-^(y)) = N, 
1=1 

where N is the number of particles (electrons)- which is considered within 
this paper as given and fixed, see (jSJ). The effective Kohn-Sham potential V 
depends on the particle densities and splits up in the following way 

V = V(u) = AE + V xc (u)-q<p, 

where (p denotes the electrostatic potential. The given external potential AE 
represents the band-edge offsets of the nanostructure materials. V xc is the 
exchange-correlation potential, which depends on the particle density, see 
Section 1. 

It is a widely discussed question how to supplement the Schrodinger op- 
erators dHD by suitable boundary conditions (see [HI H21 EES])- We choose 
homogeneous Dirichlet boundary conditions 

ip(0) = and ^(1) = 

for all ij} in the domain of the Schrodinger operator Hy. They are assisted 
by the fact that if we assume a device structure which confines the particles 
(i.e. closed system), then the particle densities vanish on the boundary of Q. 

3.2 Rigorous mathematical formulation of the problem 

In this section we give a mathematical formulation of the Kohn-Sham system; 
in particular, we make precise in which spaces the corresponding operators 
act and the solutions are chosen. In view of typical applications [6], our 
mathematical model must necessarily cover semiconductor heterostructures, 
i.e., the coefficients of the Schrodinger and the Poisson operator are in general 
discontinuous. This forecloses that the domain of the Schrodinger operator 
is not lying in W 2 ' 2 what is natural elsewhere, see e.g. [TBI US]- Fortunately, 
in the one-dimensional case the W^^-calculus already leads to satisfactory 
results. Let us introduce the Poisson operator. 
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Assumption 3.1. The dielectric permittivity e is a real, non- negative func- 
tion obeying e G L°° and ^ G 

Definition 3.2. Let Assumption 13.11 be satisfied. We define the Poisson 
operator V : Wq 1 ' 2 i — ► W~ x ' 2 by 

<Vv,w>= / £(a;)f / (a;)iy / (a;)(ix, Wj' , (17) 

where here and in the sequel < •, • > denotes the dual pairing between Wq' 2 
and W- 1 ' 2 . 

One easily estimates 

\<Vv,w>\ < \\e\\L°°\\v\\w 1 ^\\w\\w 1 ' 2 (18) 

for u,w G Wo' 2 - Consequently, P : W^' 2 h- ► W~ 1,2 is well defined and 
continuous. Furthermore, we have 

11011^.2 < ||2/e|Uo= I' e(x)\<P'(x)\ 2 dx, G W^ 2 , 
Jo 

what implies 

11011^,2 < ||2/e|U» <P0,0> 

1 2 

for G W ' . Hence, by the Lax-Milgram lemma, the inverse operator 
V~ x : W~ 1,2 i — > Wq' 2 exists and its norm does not exceed ||2/e|| £ °o. 

Assumption 3.3. The density of ionized dopants D is a 'real distribution' 
from W~ 1,2 , what means that it takes real values if applied to real elements 
from Wq' 2 . 

Definition 3.4. Let Assumptions IJ.il and \3.3\ be satisfied. Further, suppose 
u G L 1 <^-> W~ 1,2 . We say that tp G W 1,2 is a solution of the Poisson 
equation (fUj) obeying the inhomogeneous Dirichlet boundary conditions ( 173]) 
if := <p — ip G Wq ' satisfies 

V(f) = D — qu. 

We set 

^>(u):=^ + V- 1 {D-qu). (19) 

Next we are going to introduce the Schrodinger operator. 

Assumption 3.5. The effective mass m is a real, non-negative function 
obeying m G L°° and — G L°°. 
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Definition 3.6. Let Assumption 13.51 be satisfied. If V G L 1 is real valued, 
then the Schrodinger operator H v : W ' 2 i — ► W~ 1,2 corresponding to the 
potential is defined by 

h 2 r 1 1 _ r 1 _ 

< H v v,w >= — / — — v'(x) w'(x)dx + / V(x)v(x)w(x)dx, 
2 Jo rn(x) Jo 

v,w G Wo' 2 - The definition is justified, because W Q ' 2 continuously embeds 
into L°°(0, 1). Thus, the second term on the right hand side is always finite 
and defines a continuous sesquilinear form on Wq' 2 . The operator with zero 
potential will be denoted by H in the sequel. The restriction of the operators 
just introduced to other range spaces, in particular, to L 2 (0, 1), we also denote 
by H v . 

For any real valued V G L 1 (0, 1), the restriction of H v to the range space 
L 2 (0, 1) is selfadjoint and has a complete orthonormal system of eigenfunc- 
tions. All eigenvalues are then real and simple. 

Definition 3.7. We say that a continuous function / : R i— > [0, oof is from 
the class T>, if one of the following conditions is satisfied: 

(i) there is a t G] — oo, oof such that / is strictly monotonously decreasing 
on the interval ] — oo, t[ and zero on [t, oof, 

(ii) the function / obeys f(s) > for any set and 

sup f(s)s 2 < oo. 
se[i,oc[ 

is valid. 

Assumption 3.8. In the sequel we assume that all occurring distribution 
functions are from the space T>. 

Remark 3.9. Let us explicitely note that in contrast to other papers besides 
continuity no further regularity assumptions are imposed on the distribution 
functions. Only this allows to include the zero temperature case, see Chap- 
ters 6. 

Lemma 3.10. Let Assumptions VJ.5 \ and \3.8\ be satisfied. If {Xi(V)}i are the 
eigenvalues of the Schrodinger operator Hy with real potential V G L l , then 
for every N G [1, oof there is exactly one real number E G M which satisfies 

2j2f(MV)~E) = N, 

i 

This real number is called the chemical potential and is denoted by fif(V). 
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Proof. For every / eD and every E G R the expression f(Xi(V) — i?) is 
finite, lim £\ /(Aj(y)-£) = and lim £\ /;(A(\/) - E) = oo. Moreover, 

£?t— >— oo Eh- >oo 

it is not hard to see that on the set {E G R : J2i f(H v ) ~ E ) + °) the 
function E — ► ^ z /(Az(V) — -E) is strictly increasing. □ 

Remark 3.11. In the following we assume that the number of particles in 
the Kohn-Sham system is always fixed by iV G [1, oof without indicating this 
explicitly. 

Definition 3.12. Let Assumptions 13.51 and 13.81 be satisfied. Further, let 
{ipi{V)}i and {\i(V)}[ be the eigenf unctions and eigenvalues of the Schro- 
dinger operator Hy with real potential V G L 1 and let fif(V) be the chemical 
potential. Then 

M f (V)(x) :=2 J £f(\ l (V) - H (V))\MV)(x)\ 2 , xG [0,1] 
i 

defines an operator Aff : L 1 i — > L 1 which is called the particle density 
operator. 

Remark 3.13. AffiV) obviously satisfies 

-1 oo 

/ M f (V)dx = 2j2f(MV)-li f (V))=N. 
Jo i=i 

Assumption 3.14. a) The potential AE is a real-valued L l function, 
b) The exchange-correlation term in its dependence on the particle densities, 
i.e. the mapping u \ — > V xc (u), is a continuous and bounded mapping from L 1 
into L 1 . This assumption covers the Hartree-Fock-Slater and Thomas- Fermi 
exchange-correlation terms. 

Definition 3.15. Let Assumptions 13.11 13.31 13.51 13.81 and 13.141 be satisfied. 
The pair {ip, u} G W 1,2 x L l is called a solution of the Kohn-Sham system 
with distribution function /, if (p solves the Poisson equation ( JT2l) with in- 
homogeneous Dirichlet boundary conditions (|T3l) and the particle density u 
is given by 

u :=N f (AE + V xc {u) -qcp). 

4 Existence of solutions 

In this section we are going to show that the Kohn-Sham system always 
admits a solution. As in (9j [10] Schauder's fixed point theorem is used, what 
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requires several estimates (e.g. eigenvalues of the Schrddinger operator, W 1,2 - 
bounds of the particle density operator). To assure its applicability, we first 
establish some prerequisites. 



Definition 4.1. For m G L°° we set m := max ( 1 



2 m 



Remark 4.2. Recognizing that m has been defined such that 1/m is a mono- 
tonicity constant of the operator Hi = H + 1 : Wq' 2 h- > W~ 1,2 , the Lax- 
Milgram lemma shows, that the norm of the inverse operator is not larger 
than m, 



Li,2 <m< (H + l)^,^ >> IK-Ho + l) 1 1| 



By some calculations one finds 



ll L 2 , 



1.2 



which proves the continuous embedding W ( 

The following proposition allows to compare the eigenvalues of the Schro- 
dinger operators Hy and H and, additionally, provides a comparison be- 
tween the operators (Ho + l) -1 / 2 and (Hy — p)^ 1 ^ 2 - Both we will need later 
on as technical instruments. 

Proposition 4.3. Let Assumption Iff. 51 be satisfied and let V G L 1 be real- 
valued. Then the following holds: 

(i) The eigenvalues Xi(V) of the operator Hy can be estimated as follows: 

^(q + 1) + p v < A,(V) < ^(q + 1)- Pv -2, 1 = 1,2,... (20) 
where the q are the eigenvalues of the operator Hq and py is given by 

Pv = — 2||V"|||im — 1 . 

(ii) For p < py the spectrum of (Hy — p)~ l is contained in [0, 2] and 
H v -p)-%{H Q + l)% = (H + 1)Hh v - <V2. (21) 



A proof of this is to be found in [5], see Prop. 3.3. 
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Remark 4.4. The form which defines H may be estimated as follows: 
— essinf I — 1 / \v'(x)\ 2 dx < 



2 xg(o,i) [m j Jo 

— f m(x)~ 1 \v'(x)\ 2 dx < — ess sup < — 1 / \v'(x)\ 2 dx. 
2 Jo 2 xe(o,i) l m J Jo 

Thus, the eigenvalues of Ho compare by the minimax principle from below 
and above in an obvious manner with the eigenvalues of — -r^ - combined 
with a homogeneous Dirichlet condition. 

The reader should notice that Proposition 14.31 gives uniform bounds with 
respect to L 1 — bounded sets of potentials. 

From Proposition 14.31 we can deduce the following 

Lemma 4.5. Let Assumption \3.5\ be satisfied and let M. C L 1 be a bounded 
set of real-valued potentials. If 

P < Pm '■= inf pv = —2m sup ||V||xi — 1? 
then the mapping 

L l DM3V^(H v - p)- 1 e B(L 2 ) 
is Lipschitz continuous in V with a Lipschitz constant depending on A4. 

Proof. If p < pm, then p belongs to the resolvent set of Hy for any V E A4 
by Proposition 14.31 Moreover, one has 

H v -p = H v -p v + p v -p>h (22) 

since Hy — py > |J for V G A4 by Proposition 14.31 Applying the resolvent 
equation 

(Hy -p)- (Hjj -p) = (Hy -p)(U- V)(H V - p), 

one obtains 

\\(Hy-p)- l -(H U -p)- l \\ = 

WiHy-p^aj-v^Hu-p)-^. 
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The latter term may be estimated as follows: 

\\(H V - p)~\U - V)(H V - p)- 1 !! < ||E7 - V\\ BiL oo. Ll) x 
ll(^-p)- 1/2 ||||(^y-p)- 1/2 (^o + l) 1/2 ||x 

\\(H + l)- 1/2 || B (Li;L 2 )||(^0 + ir V2 |U(L-L~) X 

||(i?o + l) 1/2 (^-p)- 1/2 ||||(^-p)- 1/2 ||. 

The factors \\(Hy — p) _1 ^ 2 || and \\(Hjj — p) _1 ^ 2 || are uniformly bounded in 
U, V due to (22)). \\(H V - P y 1/2 (H + \\(H + l) l l\Hu - p)" 1 /^ are 

uniformly bounded by ( l21i) . Furthermore, \\(Hq + 1)~ 1 ^ 2 ||b(l 2 ; l° c ) is finite by 
the embedding W 1,2 <^-> This implies \\(H + 1)~ 1/2 \\b{V;I?) < oo by 

duality. Finally, \\U — V||jg(2,<x>.£i) is identical with \\U — V||x,i. □ 

Corollary 4.6. Let Assumption \3.5\ be satisfied. If {V n } n converges to V in 
L 1 , then the operator sequence {H + V n } n ^ converges in the norm resolvent 
sense to H + V . In particular, the eigenvalues and the eigenprojections of 
H + V n converge to the corresponding eigenvalues and eigenprojections of 
H + V. 

The proof follows from the preceding corollary and a well known pertur- 
bation theorem, see [TT1 Ch. IV. 3. 4]. 

Lemma 4.7. Let Assumptions \3. 5\ and \3.8\ be satisfied. 

(i) For any bounded set of real potentials in L l the set of chemical potentials 
is also bounded. Additionally, this bound can be taken even uniform 
with respect to any subset T>i C T> of distribution functions f obeying 
in addition 

sup sup f(t)t < oo and inf f\a) > for one a G R. (23) 
/ex»i te[i,oo[ / er, i 

(ii) Let C := {fj}jLi be a sequence of functions from T>\ which converges 
uniformly on bounded intervals to a function f G T> as j — > oo. If 
Vj I— >■ V in L 1 , then lim PfAVj) = Pf(V). 

j->co 

Proof, (i) By the monotonicity of /, one has J2i f(^i(V) — E) > k inf / e c f(a) 
if k items Xi — E are situated below a. This means f(\i(V) — E) i— > oo 
for E \— > oo uniform within C. Thus, the chemical potential p/(V) have to 
be bounded from above by the monotonicity of / and ( |20l) . uniformly within 
the class C. On the other hand, if E < inf 2 inf spec(Hv,) — 1, then one can 
estimate 

T f(HV) - E) < sup sup f(t)tJ2(MV)-E)-\ 

, fee te[i,oo[ , 
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what -again by (|2"U]1 - tends to zero for E i— » — oo uniformly for f E C and 
uniform over sets of potentials which are bounded in L 1 . 

(ii) First, the chemical potentials are uniformly bounded, due to (i). 
Thus, the eigenvalues of the operators ify^Av (v,) admit uniform bounds 
as in Proposition 14.31 Assume that the assertion was not true; then for a 
subsequence {Vk\k one has 

\» fk (V k ) - £ f (V)\ > e > 0. (24) 

Because the chemical potentials form a bounded set, we may again pass to 
a subsequence {Ki}n, and suppose 

\im f i fn (V n ) = E^ f i f (V). (25) 

n— >oo 

We will now show that this leads to 

lirn UHVn) - a(K)) = f^ V ) ~ E )' ( 26 ) 

I I 

First , it follows from ( 1231) and Proposition 14.31 that for any 5 > there is a 
number Iq such that 

5^/n(AKK)-/i/„(K))<5 and £/(A,(V) -£)<<$ (27) 

uniformly for all n. The remaining eigenvalues A;(K l ) for / < Iq and all 
n lie in a bounded interval and, additionally, due to Corollary 14.61 one has 
lim \i{V n ) = Xi(V) for every I. Thus, according to (|25l) . (|27|) and the uniform 

rv— >oo 

convergence of to / on bounded intervals as j ; — > oo, the term 

E/nC^TO -maCK)) -^)l 

becomes smaller than 35 for sufficiently large n and arbitrarily chosen 5 > 0. 
This implies (1261) . but (1261) must be false: the terms on the left hand side 
all equal N, what cannot be true for the right hand side due to (j23|) and 
Lemma [3.101 Hence, (j2*4l is wrong, what proves (ii). □ 

Remark 4.8. The lemma shows in particular that the chemical potentials 
continuously depend on the potential. 

Theorem 4.9. Let Assumptions \3.5\ and \3.8\ be satisfied and let M. be a 

bounded set of real potentials in L 1 . Then the following holds: 
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(i) The image Aff (Ai), f GD, is a bounded set in W ' 2 . The bound may be 
taken uniformly with respect to any set V 2 C V of distribution functions 
f which satisfy the additional conditions 

sup sup f(t)t(t + 1) < oo and inf f(a) > for one a6l 

/ex> 2 te[o,oo[ / eX, 2 

(28) 

(ii) The particle density operator Aff : L 1 \— > L 1 , f <EV, is continuous. 
Proof, (i) For V € M we get 

||^/(V0|| w i., = ||2^/(A I (lO-A*/(V r ))|^ll^< (29) 
i 

2 f(Hv) - ii f {v))\\ WWwf ^ 4 E - 



where in the last step we used |||'0| 2 |lw rl ' a - 2 llV ; llLi,2 for all ?/> £ W^ 1,2 . We 
estimate the terms in (1291): 



IWI^<mll(flo + Wi||^< 

m||(# + l)*(flv - p)"* flK^v - < 2m(\ l (V) - p), 

(see ( 121]) ) where p < p^ where pm given by Lemma 14.51 is a uniform lower 
bound for the spectra of the operators Hy with V £ A4. Thus, ( 1291) is not 
larger than 

Wf(V)\\ w w < 8mJ2f(MV) - t*f{V))(\i(V) -p)< 
i 

8m £ /(W) ~ m/(*0) - p) + (W) - p/(v))] < 

8miV| M/ (V-) - p| + 8m £ /(A,(V) - M/OOXW) - p f (V)). 

l:\i(V)-n f (V)>0 

The last sum may be estimated by 

sup f(t)t(t + l) Yl (M^) - + I)" 1 . 

' e[0 '°° [ *:AiOO-/ff( v ')> 



Obviously, the condition ( 1281) is stronger than ( 1231) ; thus, the set of chemical 
potentials is uniformly bounded for / £ C and V £ A4. This, together with 
the eigenvalue estimates (12"U1) and Remark |4~41 proves (i). 
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(ii) According to Corollary 14.61 and Lemma 14.71 the eigenvalues and the 
corresponding chemical potentials depend continously on V 3 L 1 . Further- 
more, if {Vj}j converges in L 1 to V, then - due to Corollary 14.61 - the (one- 
dimensional) eigenprojections Pj^ = (• , fj,k) fj,k converge to the correspond- 
ing eigenprojections Pj~ = {-,fk)fk- Applying these eigenprojections to the 
vector f k , one easily obtains 

\(fk,fj,k)\ ^ 1 and (f k Jj,k)fj,k^fk inL 2 for j i-> oo. 

From this it is not hard to see that \fj t k\ 2 l— * \fk\ 2 m L 1 . Observing that for 
sufficiently large M 

II f(W) - nf(v)) \HV)\ 2 \W < £ f^M - vAv)) 

1>M 1>M 

can be made arbitrarily small uniformly over an L 1 -bounded set of potentials 
V, one obtains (ii). □ 

Corollary 4.10. Let Assumptions \3.5\ and \3.8\ be satisfied. If X is a Banach 
space continuously injecting into L 1 such that W ' 2 compactly embeds into 
X , then the particle density operator A/} : L 1 i— > X is well defined and 
continuous. 

Proof. That it is well defined follows from the embedding Wq' 2 X. As- 
sume that the continuity does not hold; then there must be a sequence {V n } n 
converging in L 1 to V such that 

||A/}(K)-A/}(V)||x>e>0. (30) 

The statement (i) of the foregoing lemma tells us that {A//(V r ri )} n is bounded 
in W 1,2 . Thus, by the compactness of the embedding W 1,2 X there must 
be a subsequence {Vk]k such that A//(Vfc) converges in X to an element W. 
But A//(Vfc) converges to A//(V) in L 1 . Thus W must equal Mf(V) by the 
continuous injection X L l , what contradicts (1301) . □ 

To prove existence of solutions of the Kohn-Sham system we will introduce 
an appropriate subset of L 1 together with a suitable mapping $ from this set 
into itself. $ will be constructed such that the solutions to the Kohn-Sham 
system coincide with the fixed points of this mapping. 

Definition 4.11. We set 

L l N := I u E L l : u > 0, / u dx = N 



o 
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Let Assumptions 13.11 13. 3\ \3.5\ 13.81 and 13.141 be satisfied. Then we define the 
mapping ® f : L l N -> L l N by 

:= Aff(AE + V xc (u) - q<p(u)) (31) 

where <£>(tt) denotes the w-dependent solution to Poisson's equation given by 

una. 



The task now is to verify that has a fixed point. We will use Schauder's 
theorem to achieve this. 

Theorem 4.12. Let Assumptions \3J\ [OJ GL21 and \3J4\ be satisfied. 
Then 

(i) the pair {(p(u), u] is a solution of the Kohn-Sham system if and only if 
u G L l N is a fixed point of <3>/ ; where ip(u) is given by / TTPj) 

and 

(ii) the mapping $ f has a fixed point. 

Proof, (i) The first part of theorem follows immediately from the definition 
of the mappings A/} and 

(ii) To prove the second part we note that L l N is a closed, bounded and 
convex set, which is mapped by into itself (Def. 14.111) . 

Continuity: Since L 1 <— ► W~ 1,2 , the solution (p to Poisson's equation 
depends continuously (in W ' 2 ) on u. Hence, due to Assumption 13.141 the 
mapping L 1 3 u \— > AE + V xc (u) — qf{u) G L 1 is continuous. Theorem 14.91 
then implies the continuity of 

Compactness: According to Theorem I4.9[ the image of a L 1 -bounded 
set of potentials is bounded in the space W ' 2 . The compactness of the 
embedding W 1,2 L 1 yields the compactness of $/. 

Thus, due to Schauder's fixed point theorem, has a fixed point in 
L l N . □ 



5 Particle density operator and monotonicity 

In this section we want to show some additional properties of the particle 
density operator. For this it is necessary to restrict some considerations to 
the real parts of spaces which were up to now considered as complex ones. We 
indicate this by an additional subscript R, e.g. L§?. The upcoming results 
are based on the following theorem. 
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Theorem 5.1. Let H be a self- adjoint operator in the separable Hilbert space 
Tt with compact resolvent and let U and V be bounded, self-adjoint operators 
onTi. If f : R — > R is a Borel measurable function such that f(H + U) and 
f(H + V) are trace class operators, then the formula 

tx([f(H + U)- f(H + V)](U - V)) = (32) 

oo 

X)(/(A*)-/0*i))(Aik-Mi)K^,6>l 2 

k,l=l 

is valid. Here {Xk}k is the sequence of eigenvalues of H + U (H + 

V) and {ipk}k ({£z}j) i s an ~~ orthornormalized - sequence of corresponding 
eigenvectors. 

The proof is given in [8]. 

Corollary 5.2. If Assumptions fX31 and \3^R are satisfied, then the mappings 
L R 3 V -> A//(V) E and Wjj 3 V —> A/}(V) G W^ 1 ' 2 are monotone. 

Proof. We specify the Hilbert space H to L 2 and first consider potentials 
U, V e Lj^P, which are identified with the induced multiplication operators 
on L 2 . Replacing U, V in the preceding theorem by U — Hf{U), V — Hf{V) 
and observing 

tr([f(H + U- N {U)) -f(H + V- iMf{V))](n f (U) - H {V))) = 
tr([f(H + U- H {U)) -f(H + V- iif{V))\){^{U) - p f (V)) = 

one obtains 

J M(U)-tf,(V))(U-V)dx = 

2tx([f(H + U- fM,(U)) -f(H + V- n f (Y))](U - V)) = 

2tr([/(ff + U- n f (U)) -f(H + V- /i f (V))](U - H (U) - V + ^ f (V))). 

In view of fl32|) the right hand side is negative because the distribution func- 
tion / is monotonously decreasing. 

Let now U, V e L\ be arbitrary. If {U n } n , {V n } n are sequences from L^ 
which converge to U, V in L 1 , respectively, then, on one hand, Aff{U n ) i— > 
M f (U) and Mf{y n ) i-> Af f (V) in Lf, due to Corollary Offl Because one 
already knows that /„ '(Af f (U n ) - N f (V n )){U n - V n ) dx < 0, one obtains 
tiWf{U)-M f {V)){U -V)dx<0. 

2 12 i 

The W 0R <-y W R ' duality extends the L R <-> L R duality; thus, the 
second assertion follows from the first. □ 
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Theorem 5.3. Let Assumptions\3J\\EM IX51EQ1 and \3.14\ a) be satisfied. If 
the correlation term V xc is absent, then the Kohn-Sham system has a unique 
solution {ip, u}. 

Proof. It is not hard to see that in this case the system can be written as 
one equation for the electrostatic potential in the real space W R ' , namely 

- ^~7^<t> + gA/}(A£ -q$-q<P) = D (33) 

where ip = <p + 0. Since the operator — '■ W Q ^ i— > W-^ ' is strongly 

monotone, the operator —^-^ + qJ\ff(AE — qip — q ■) : \— > ' 2 is 
also strongly monotone by the foregoing corollary. Additionally, this latter 
operator is continuous, and hence (1331) has a unique solution by the theory 
of monotone operators, see [21] Ch.26.2. □ 

Corollary 5.4. Let the assumptions of Theorem \5.3\ be satisfied. If{<p,u} is 
a solution of the Kohn-Sham system, then the electrostatic potential = (p—(p 
satisfies the following a priori estimate: 



w w < ■^(||D|| w --i,2 + 7 il . w -i, a JVg). 



11 12 

where M is the monotonicity constant for the operator — '■ Wok ^ 

WjgT 1,2 and 7_li ; vk- 1 > 2 is the norm of the embedding operator L 1 W^ 1 ' 2 . 
Proof. Clearly, satisfies the equation 

-^£^0 + qM f {AE + V xc {u) -qtp- q<f>) = D, 

which can be regarded as an equation in W R ' , due to Assumption 13.31 
and the fact that the dielectric permittivity matrix e has real entries, cf. 
Assumption 13.11 Considering u as fixed, the left hand side is a strongly 
monotone, continuous operator when acting on 0. We denote it for brevity 
by A. Its monotonicity constant is at least M. Using the strong monotonicity, 
we may estimate: 

\\M-M\wM\<i>\\w^ > I < -40 - -40, > | > M||0||^, 2 , 

M 0,K KK 0,IR 

what leads to 

W\\ w g < ^(\\D\\ w -^ + q\\M f (AE + V xc (u)-qm w ^) < 
< Jj(\\ D Ww^ + W L i;W^\\M f (AE + V xc (u) - q$)\\ L i 

Obviously, \\Aff(AE+V xc (u) — q(p) H^i equals N, what gives the assertion. □ 

Remark 5.5. Note that this estimate depends in no way on the distribution 
function / (within the class of monotonously decreasing functions, of course). 
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6 Convergence to zero temperature 

In the following we introduce the function 9 by defining 

0{x) = I 1 : X ~ ] 
[ x : x > 1. 

We start with the following technical lemma. 

Lemma 6.1. Let Assumptions Iff. 51 and \3.8\ be satisfied. Further, let C : = 
{fj}fLi c ^2? cf. Theorem \4-9[ and f G V such that 

lim sup \fj(x) -f(x)\9(x) = (34) 

■J 1 -* 00 x£[a,oo[ 

holds for every a g] — oo, — 1[. If {Vj}°? =1 , Vj G L 1 converges to the real 
potential V G L 1 in L 1 , then lim A/V (V,-) = Aff(V) in L 1 . 

Proof. One has the estimate 

IM,M)-WS)IU* = 

sup | l{M fj {V j )-M f {V j ))Wdx\ = 

\\W\\ L °°=l J 

sup |2tr ((fiiH Vi - /^(V})) - f(H Vj - p f (Vj))) W) \ < 

\\W\\ L oc=l 

2\f 3 (H Vj -HfM)) -f(H v . - MVS-)) Ik < 

2 l(/j(^v r j - Vfj(Vj)) - f(Hvj - V>fj(Yj)))(H Vj - p)\\g\\(H Vj - ^) _1 ||bi, 
where p < pc, cf. Lemma [4.51 This leads to the estimate 
ll^(^-A/}(^|| L i< 

Wi(H Vj -nM)) - f(H Vj - PfMMHvi - P )h\\{H Vj - p)- 1/2 \\ 2 B2 , 
Using (1211 . one estimates the second factor by 

II(^-p)- 1/2 IIb 2 < 

\\(H Vj - p)- l ' 2 (H - p) 1/2 h\\(H - p)-y*\\ B2 < V2\\(H - P)- 1/2 W 
To estimate the first one, we write 

\\{f 3 {H V] - p fj (V 3 )) - f(H Vj - p f (Vj)))(H Vj - P )\\b< 



sup 

t£[infj spec(Hyj),oo[ 






-f(t 




-p)< 




sup 

tE[infj spec(Hyj),oo\ 


\(fj(t 




-in 




-p) + 


(35) 


sup 

tE[infj spec(Hyj),oo\ 


\(f(t- 


-H 3 m- 


-f(t- 




p)- 


(36) 
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The term (|35j) converges to zero due to flMl) and CCP 2 . 

We consider the term flHUj) . Lemma 14.71 yields that {^f J (Vj)}'^ =1 and 
{/i/(V 3 )}°? =1 converge to /z/(V). Thus, the restriction of /(• (Vj)) — /(• — 
/i/(Vy)) to finite intervals uniformly converges to zero by the continuity of 
/. On the other hand, / decays at oo as ^, thus, for large arguments t the 
absolute value of f(t — fif.(Vj)) — f(t — fif(Vj)) becomes arbitrarily small 
uniformly in j. This altogether shows that (1361) goes to zero. □ 



Corollary 6.2. Let Assumptions [3J\ [X3| \EM E3 and \3.14\ be satisfied. 



Further, let C := {fj}'jL 1 ^ ^2; c f- Theorem \4-9j and f G T> such that 
condition (5^] ) holds for every o g] -oo,-1[. If Uj G Ljy converges 

to the real function u G L X N in L 1 , then lim $/■(«,) = in L l . 

Proof. By definition fl3~T|) we have 

$ /. ( Uj ) := JV^ (AE + V xc {uj) - ^(a,)), «j G Ljv, j = 1, 2, .... , 
where 

<^(t%) = & + V^ 1 (D - quj). 
Since lim ip(iij) = (p(u) = <p+V~ l (D—qu) in W 1 ' 2 one gets that lim <p(itj) = 
(p(u) in L l . Therefore lim V, ■ = V in L 1 , where 

Vj := AE + V xc (uj) - q<p(uj) and V := AE + V xc {u) - qip(u). 



Applying Lemma [6. 1[ we complete the proof. □ 
Theorem 6.3. Let Assumptions \3J\ [P. [Q. LOI and 6e satisfied. 



Further, let f G £> and C := {/j}^ ^ P 2; c/. Theorem \4-9\ such that 
C := {fj} ( jL 1 obeys (34\) for every a g]— oo, — 1[. If{{(pj, Uj}}JL 1 are solutions 
of the Kohn-Sham system with respect to the distribution function fj, then 
there is a subsequence {{<fk, u k}}'k=i which converges in L°° xL 1 to a solution 
{ip,u} of the Kohn-Sham system with distribution function f . 

Proof. By Theorem 14. 121 {ipj , Uj} G W 1 ' 2 x L l is a solution of the Kohn-Sham 
system with respect to the distribution function fj if and only if Uj = (uj), 
j = 1,2,... and the corresponding potential is given by <pj = ip(v,j) = <p + 
V~ 1 (D — quj). According to Theorem 14.91 there are subsequences {m*;}^ 
and {ip^kLi suc h that the following properties are satisfied: 



The sequence {wfc}£° is bounded in W 1 ' 2 , obeys Uk G L l N and converges 

-l 



in L 1 to an element u G L) 
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• The sequence of potentials {<Pk}kLi converges in L 1 and, additionally, 
weakly in W 1,2 to an element (p. 

By Theorem 14.121 the pair {ipk,Uk} is a solution of the Kohn-Sham system 
with distribution function ff. if and only if Uk is a fixed point of the map 
i.e. 

u k = ®f k {u k ), k = 1,2, . . . . 
and the potential <pk is given by 

¥k = ^f> + V' l {D - qu k ). 

By lim Uk = u in L 1 and Corollary 16.21 we get u = for u & L l . This 

shows that w is a fixed point of Moreover, one has <p = ip + lim V~ X (D — 

quk) in L°° which shows that <p G ly 1 ' 2 . By Theorem 14.121 the pair {(p, u} is 
a solution of the Kohn-Sham system with distribution function /. □ 

If the Kohn-Sham system with distribution function / has several solu- 
tions, then it remains unclear to which of them a sequence of solutions of 
Kohn-Sham systems with distributions functions fj converges. However, if 
the exchange correlation term is absent, then the result can be improved. 

Corollary 6.4. Let the assumptions of Theorem \6.^\ be satisfied. If the ex- 
change correlation term V xc is absent and if {{y^j; are unique solu- 
tions of Kohn-Sham systems with distribution function fj, then {{<Pj, Uj}}^° =1 
converges in L°° x L 1 to the unique solution {ip, u} of the Kohn-Sham system 
with distribution function f . 

Proof. Assume that the sequence {{<Pj, Uj}}'jL 1 does not converge to {ip,u}. 
In this case there is a subsequence {{<Pk,Uk}}kLi converging in L°° x L 1 to 
an element {<p, u] G L°° x L 1 which is different from {ip, u}. However, by 
Theorem 16.31 the pair {p, u] is a solution of the Kohn-Sham system with dis- 
tribution function /. Since this Kohn-Sham system admits only one solution 
the solution {<p, u} coincides with {ip, u} - what is a contradiction. □ 

Lemma 6.5. The function fp{x) = 4 ln(l + e~^ x ), x G K ; strictly decreases 
in (3 G]0, oof. 

Proof. One calculates 

f) I fv* p 
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what immediately shows the assertion for x > 0. Putting — fix =: 7, the 
assertion for negative x is equivalent to 

■P^— < ln(l + e 7 ) , 7 > 0, 

which follows from < 7 = ln(e 7 ) < ln(l + e 7 ). □ 

In order to apply this to the Kohn-Sham system at zero temperature, we 
show in the following that the corresponding distribution function satisfies 
the condition (SI 



Lemma 6.6. Let {7}}°^ be any positive sequence converging to zero. We 
set fj(x) := J-hi(l + e _/3jX ) where ftj = ^r. Further, we set 

'<*> == { 1 ; « > 0. < 37 > 

TTien condition (34\ ) is satisfied. 
Proof. We have 

lim sup \fj(x) — f(x)\0(x) < 

lim sup \fj(x) — f(x)\ + Hm sup \fj(x) — f(x)\x (38) 

J Woo a<s<l i HO °i>l 

for a < —1. Obviously we have 



-ln(l + e -^)x<^— x>l. 



This shows that the second term of (158]) tends to zero as j ' ^ 00. Further, it 
is almost obvious that fj(x) converges pointwise to the continuous function 
/ for j 1— > 00. Because the family {fp}p is monotonously decreasing in ft 
by the preceding lemma, the convergence is uniform on bounded intervals 
by Dini's theorem. This proves that the first term of (|38|) tends to zero as 
j — > 00. □ 



Acknowledgment 

H.C, H.N. and J.R. acknowledge support from the Danish F.N.U. grant 
Mathematical Physics and Partial Differential Equations. This work was 
initiated during a visit of K.H, H.N. and J.R. at the Department of Math- 
ematical Sciences of the Aalborg University and they are thankful for kind 
hospitality extended to them during the work on this paper. 



24 



References 



[1] Olivier Bokanowski, Jose L. Lopez, Oscar Sanchez, and Juan Soler. Long 
time behaviour to the Schr6dinger-Poisson-X a systems. In Mathematical 
physics of quantum mechanics, volume 690 of Lecture Notes in Phys., 
pages 217-232. Springer, Berlin, 2006. 

[2] Olivier Bokanowski, Jose L. Lopez, and Juan Soler. On an exchange in- 
teraction model for quantum transport: the Schrodinger-Poisson-Slater 
system. Math. Models Methods Appl. Sci., 13(10):1397-1412, 2003. 

[3] M Cahay, M. McLennan, S. Datta, and M. S. Lundstrom. Importance 
of space-charge effects in resonant tunneling devices. Appl. Phys. Lett., 
50(10):612-614, 1987. 

[4] P. A. M. Dirac. Note on exchange phenomena in the Thomas atom. 
Proceedings Cambridge, 26:376-385, 1930. 

[5] W. R. Frensley. Boundary conditions for open quantum systems driven 
far from equilibrium. Rev. Modern Phys., 62:745-791, 1990. 

[6] H. Gajewski. Analysis und Numerik von Ladungstransport in Halbleit- 
ern. Report 6, WIAS, Berlin, 1993. 

[7] P. Hohenberg and W. Kohn. Inhomogeneous electron gas. Phys. Rev., 
136(3B):B864-B871, 1964. 

[8] H.-Chr. Kaiser, H. Neidhardt, and J. Rehberg. Monotonicity proper- 
ties of the quantum mechanical particle density. Preprint 1275, WIAS, 
Berlin, 2007. 

[9] H.-Chr. Kaiser and J. Rehberg. About a one-dimensional stationary 
Schrodinger-Poisson system with Kohn-Sham potential. Z. Angew. 
Math. Phys., 50(3):423-458, 1999. 

[10] H.-Chr. Kaiser and J. Rehberg. About a stationary Schrodinger- 
Poisson system with Kohn-Sham potential in a bounded two- or 
three-dimensional domain. Nonlinear Anal., 41(1-2, Ser. A: Theory 
Methods) :33-72, 2000. 

[11] T. Kato. Perturbation theory for linear operators. Die Grundlehren der 
mathematischen Wissenschaften, Band 132. Springer- Verlag New York, 
Inc., New York, 1966. 



25 



[12] T. Kerkhoven. Mathematical modeling of quantum wires in periodic 
heterojunction structures. In Semiconductors Part II, volume 59 of 
The IMA Volumes in Mathematics and its Applications, pages 237-253. 
Springer, New York, 1994. 

[13] T. Kerkhoven. Numerical nanostructure modeling. Z. Angew. Math. 
Mech., 76(Suppl. 2):297-300, 1996. 

[14] W. Kohn and L. J. Sham. Self-consistent equations including exchange 
and correlation effects. Phys. Rev., 140(4A):A1133-A1138, 1965. 

[15] Elliott H. Lieb and Barry Simon. The Hartree-Fock theory for Coulomb 
systems. Comm. Math. Phys., 53 (3): 185-194, 1977. 

[16] P.-L. Lions. Solutions of Hartree-Fock equations for Coulomb systems. 
Comm. Math. Phys., 109(1) :33-97, 1987. 

[17] N. D. Mermin. Thermal properties of the inhomogeneous electron gas. 
Phys. Rev., 137(5A):A1441-A1443, 1965. 

[18] F. Nier. Etude mathematique et numerique de modeles cinetiques quan- 
tiques issus de la physique des semi-conducteurs. Master's thesis, Ecole 
Polytechnique, 1991. 

[19] F. Nier. A variational formulation of Schrodinger-Poisson systems in 
dimension d < 3. Comm. Partial Differential Equations, 18(7-8):1125- 
1147, 1993. 

[20] E. Prodan and P. Nordlander. On the Kohn-Sham equations with peri- 
odic background potentials. J. Statist. Phys., lll(3-4):967-992, 2003. 

[21] P. N. Racec, U. Wulf, and J. Kucera. Integration of quantum transport 
models in classical device simulators. Solid-State Electronics, 44(5) :881- 
886, 2000. 

[22] Oscar Sanchez and Juan Soler. Long-time dynamics of the Schrodinger- 
Poisson-Slater system. J. Statist. Phys., 114(1-2): 179-204, 2004. 

[23] J. C. Slater. A Simplification of the Hartree-Fock Method. Phys. Rev., 
81(3):385-390, 1951. 

[24] E. Zeidler. Nonlinear functional analysis and its applications. II/B. 
Springer- Verlag, New York, 1990. Nonlinear monotone operators, Trans- 
lated from the German by the author and Leo F. Boron. 



26 



B. Zimmermann, E. E. Marclay, M. Ilegems, and P. Gueret 
Self-consistent calcultions of tunneling currents in n + -GaAs/i 
AL r Gai_ a As/n + -GaAs structures and comparison with measurements 
J. Appl. Phys., 64(7):3581-3588, 1988. 



27 



